Stable Vortex Solitons in Nonlocal Self-Focusing Nonlinear Media 
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We reveal that spatially localized vortex solitons become stable in self-focusing nonlinear media 
when the vortex symmetry-breaking azimuthal instability is eliminated by a nonlocal nonlinear 
response. We study the main properties of different types of vortex beams and discuss the physical 
mechanism of the vortex stabilization in spatially nonlocal nonlinear media. 
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Vortices arc fundamental objects which appear in many branches of physics 0. In optics, vortices are usually 
associated with phase singularities of diffracting optical beams, and they can be generated experimentally in different 
types of linear and nonlinear media However, optical vortices become highly unstable in self- focusing nonlinear 
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, media due to the symmetry-breaking azimuthal instability, and they decay into several fundamental solitons Q- In 
spite of many theoretical ideas to stabilize optical vortices in specific nonlinear media no stable optical vortices 
created by coherent light were readily observed in experiment [3. Thus, the important challenge remains to reveal 
physical mechanisms which would allow the first experimental observation of stable vortices in realistic nonlinear 
media. 

In this Letter, we reveal that the symmetry-breaking azimuthal instability of the vortex beams can be eliminated 
in a medium characterized by a nonlocal nonlinear response. This observation allows us to suggest a simple and 
^ , realistic way to generate experimentally the first stable spatially localized vortices in self-focusing nonlinear media. 
We study the main properties and stability of different types of vortex beams, and discuss the physical mechanism of 
their stabilization in spatially nonlocal nonlinear media. 

We notice that there are many physical systems characterized by nonlocal nonlinear response. In particular, a 
C^l I nonlocal response is induced by heating and ionization, and it is known to be important in plasmas [6|. Nonlocal 
^ . response is a key feature of the orientational nonlincaritics due to long-range molecular interactions in nematic liquid 
■ crystals . An interatomic interaction potential in Bose-Einstein condensates with dipole-dipole interactions is also 
substantially nonlocal Q . In all such systems, nonlocal nonlinearity can be responsible for many novel features such 
as the familiar effect of the collapse arrest |3, [l3 ■ 

We consider propagation of the electric field envelope E{X, Y, Z) described by the paraxial wave equation: 



where ko is the wave number, and the function O characterizes a nonlinear, generally nonlocal, medium response. 
For example, in the case of the wave beam propagation in partially ionized plasmas, O = T^/T is the relative 
J> . electron temperature perturbation, T being unperturbed temperature, the coupling coefficient tit = +1- Stationary 
temperature perturbation obeys the equation 
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where ~ 3mT([jQ -f v'^)/e'^, v^, is the electron collision frequency, cjq is the wave frequency, m is electron mass, 
c? characterizes the relative portion of the energy that electron deliver to a heavy particle during single collision. 
The second term describes thermal diffusion with the characteristic spatial scale . Note, that the model identical to 
Eqs. 1^ and (j^J has been employed in Ref. |3| to study two-dimensional bright solitons observed experimentally in 
nematic liquid crystals 0,^3- this case, the field 8 describes the spatial distribution of the molecular director. 

Rescaling the variables, {X,Y) = le{x,y), Z = 2lez/e and the fields, E = {E^e/ ^/rLT)'^{x,y,z) and Q — 
{t^ /nT)0{x,y, z), where e ~ (fco^e)"^; we present Eq. (jSJ in the dimcnsionless form. 
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where Aj^ = d"^ /dx^ -f- /dy^ is the transverse Laplacian. For the analysis performed below, we omit in Eq. Q the 
term proportional to . 
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Thus, the basic dimensionlcss equations describing the propagation of the electric field envelope '^{x,y, z) coupled 
to the temperature perturbation 9{x,y, z) become 

a* 

dz (4) 



a 



A, 6*= IfP. 



In the limit ^ 1, we can neglect the second term in the equation for the field 9 in Eq. I^J and reduce this system 
to the standard local nonlinear Schrodinger (NLS) equation with cubic nonlinearity. The opposite case, i.e. ^ 1, 
can be referred to as a strongly nonlocal regime of the beam propagation. 

We look for stationary solutions of the system Q in the form, ^(x,y,z) = ip^r) e'xjpiimip + lAz), where (p and 
r = y^E^ + y^ are the azimuthal angle and radial coordinate, respectively, and A is the beam propagation constant. 
Such solutions describe either the fundamental optical soliton, when to = 0, or the vortex soliton with the topological 
charge m, when m ^ 0. 

The beam radial profile smd accompanied temperature field 9{r) can be found by solving the system of ordinary 
differential equations, 

(5) 



where Ar'"'' = (P/dr^ + {\/r){d/dr) — [m? /r^), and we rescale ■0j '^/f'^, and by the factor A, and choose the 
unity propagation constant. Boundary conditions are: for the localized vortex field, ip{oo) = "0(0) = 0, and for the 
temperature field, d9/dr\r=o = and 9{oo) = 0. 

The second equation of the system (0J| can be readily solved for axially symmetric intensity distribution I'fp, 
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where Go is the Green's function defined at = from the general expression 

G.(ei,6,«) ~ I l^(^a^,)K,{a^,l 6 < Ci < +oo, 

and lu and K^, are the modified Bessel functions of the first and second kind, respectively. Thus, Eqs. (jSJ are equivalent 
to a single integro-differential equation obtained from jSJl when 6{r) is eliminated, or to a single integral equation: 
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^{r)= 0(77)V'(77)G„,(r, 77; 71)77 dr?, (8) 







where Gm is defined by Eq. (0) and 9 is given by Eq. 

We solve the nonlinear integral equation (jS)) using the stabilized relaxation procedure similar to that employed in 
Ref. Figure n shows several examples of the solutions of the system (0) found numerically for different values 

of nonlocality parameter a. To characterize these solutions, we define the effective radii and rg of the intensity 
distribution |7/;|^ and the temperature perturbation distribution 9, respectively, as follows, 

1 /■ 2, M2^2 2 / r^O{r)d''v 
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Figure EJa) shows the radii r^, and re as functions of the nonlocality parameter a. Both and rg decrease monoton- 
ically when the nonlocality parameter grows. In the local limit (a ^ 1), r^ and rg saturate at the same finite value, 
which increases with the topological charge. Figure [21[b) shows the beam power P ~ J \E\^ d^v as a function of the 
nonlocality parameter a. 

The important information on stability of the vortex solitons can be obtained from the analysis of small pertur- 
bations of the stationary states. The basic idea of such a linear stability analysis is to represent a perturbation as 
the superposition of the modes with different azimuthal symmetry. Since the perturbation is assumed to be small, 
stability of each linear mode can be studied independently. Presenting the nonstationary solution in the vicinity of 
the stationary mode as follows. 
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e(r,z) = e{r)+§+{r)e 



"'+'^'f + r_{r)e 



where |e±| ^ -0, \d±\ 6, ip, 6 are assumed to be real without loss of generality, we linearize Eqs. (0J) and obtain the 
system of linear equations of the form: 



The Hankel spectral transform was applied to reduce the integro-differential eigenvalue problem (O to linear alge- 
braic one. The maximum growth rate |Imu;| of linear perturbation modes is shown in Fig. |2Ia) for a single-charge 
vortex (m = 1). The symmetry-breaking modes can become unstable only for L = 1, 2, 3. The growth rates saturate 
in the local regime a S> 1. The largest growth rate as well as the widest instability region has the azimuthal mode 
with the number L = 2. The real and imaginary parts of the eigenvalues lu for this most dangerous mode are shown 
in Fig. I^Jb). Importantly, there exists a bifurcation point acr ~ 0.12 below which the growth rate Im (w) vanishes. 
Thus, the symmetry-breaking azimuthal instability is eliminated in a highly-nonlocal regime: all growth rates vanish 
provided a < aa- 

We also perform the linear stability analysis for multi-charge vortices with the topological charges m = 2 and 
m = 3. Figure ini^c) shows the growth rate of the linear perturbation modes for the vortex with m = 2. Importantly, 
the growth rate of the L ^ 2 mode always remains nonzero, and the same result holds for the vortices with m = 3. 
Therefore, the linear stability analysis predicts the existence of stable single-charge vortex in a highly nonlocal regime, 
while the multi-charge vortices are shown to be unstable with respect to decay into the fundamental solitons, even in 
the limit a ^ 0. 

Our stability analysis has been verified with direct simulations of the propagation dynamics of perturbed vortex 
solitons by employing the split-step Fourier method to solve Eqs. numerically. The results of our linear stability 
analysis agree well with the numerical simulations. In particular, the symmetry-breaking instabilities have been 
observed in the region predicted by the linear stability analysis, and some examples of the vortex decay instability are 
presented in Fig. 0] for the vortices with m = 1. If a perturbation is applied to a single-charge vortex in the strongly 
nonlocal regime (such that all azimuthal instabilities arc completely suppressed by the nonlocality) , the vortex beam 
evolves in a quasi-periodic fashion: the effective radii and amplitudes oscillate with z. Thus, our numerical simulations 
indicate that single-charge vortex solitons become stable if the nonlocality parameter a is below some critical value 
which is very close to the value acr ~ 0.12 predicted by the linear stability analysis. In experiments, the input beam 
may differ essentially from the exact stationary solution. Therefore, we perform additional numerical simulations for 
singular Gaussian input beams of the form ^'(r, 0) = /ir exp(— r^/w^ -I- iip) and, as follows from Fig.jS] observe that 
such beams are indeed stabilized below the critical value of nonlocality, and the beam effective intensity, defined as 
= J \'ii\'^(Pr, underdoes large-amplitude oscillations. 

The physical mechanism for suppressing the symmetry-breaking azimuthal instability of the vortex beam in a 
nonlocal nonlinear medium can be understood as being associated with effective diffusion processes introduced by a 
nonlocal response. Indeed, if a small azimuthal perturbation of the radially-symmetric vortex deforms its shape in 
some region, the corresponding temperature distribution along the vortex ring becomes nonuniform. As a result, the 
intrinsic thermodiffusion processes would smooth out this inliomogeneity and suppress its further growth, leading to 
the complete vortex stabilization in a highly-nonlocal regime. 

In conclusion, wc have studied the basic properties and stability of spatially localized vortex beams in a self- 
focusing nonlinear medium with a nonlocal response. We have found that single-charge optical vortices can be 
stabilized with respect to any symmetry-breaking azimuthal instability in the regime of strong nonlocal response, 
whereas multi-charged vortices remain unstable and decay into the fundamental solitons for any degree of nonlocality. 
Our numerical simulations confirm that the stable propagation of single-charge vortex is indeed possible in a nonlocal 
nonlinear medium. We expect that these results will stimulate the first experimental observation of stable vortices in 
self-focusing nonlinear media. 
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FIG. 1: Examples of the stationary vortex solutions for m = 1 and m = 2 at different values of the nonlocality parameter a. 
Shown are the fields '!/'(r) (solid) and 6{r) (dotted). 




FIG. 2: (a) Effective radii of the field intensity distribution (solid) and the temperature field rg (dotted) vs. the nonlocality 
parameter a, for m = 1,2. (b) Power P vs. a. 
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FIG. 3: Maximum growth rate of linear perturbation modes vs. the nonlocality parameter a for the vortices with (a) m = 1 
and (c) m = 2. Numbers near the curves stand for the azimuthal mode numbers L. (b) Real and imaginary parts of the 
eigenvalue ui of the most dangerous azimuthal mode with L — 2 and ?n = 1. 
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FIG. 4: Evolution of the beam intensity l^fp (upper row) and temperature field 6 (lower row) of a perturbed single-charge 
vortex soliton for (a.) a — 5 and (b) a = 1. 
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FIG. 5: (a) Stable propagation of a single-charged vortex generated by a singular Gaussian beam with m = 1 {h — 0.14, 
w = 4.2, a = 0.07) shown for the beam intensity l^*]^ and temperature field 9. (b) Oscillatory dynamics of the amplitude of 
the intensity field vs. z during the beam propagation. 



